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I.  INTRODUCTION 


As  shown  by  Epstein  et  al.  (1980;  hereafter  referred  to  as  Paper  I), 
single-antenna  measurements  of  the  combined  signal  from  Saturn  and  its  rings 
made  over  many  years,  as  the  ring  inclination  changes,  can  be  used  to 
determine  separate  values  of  the  ring  temperature  and  the  disk  temperature. 
However,  two  factors  limited  the  specificity  of  Paper  I's  results:  the  large 
errors  of  the  laboratory  measurements  of  the  absorption  coefficient  of  water 
ice  (a  likely  major  constituent  of  the  ring  particles),  and  the  large  error  of 
the  derived  ring  brightness  t  *verature  because  of  the  lack  of  observations  at 
small  ring  inclinations.  Th  ^aper  utilizes  new  laboratory  ice  measurements 
(Mlshima  et  al.,  1983)  and  low- inclination  observations  to  derive  more 
accurate  ring  properties. 


II.  OBSERVATIONS  AND  RESULTS 


Observations  at  low  ring  inclinations  were  made  at  90  GHz  (3.3  mm)  with 
the  Aerospace  4.6-m  radio  telescope,  as  described  in  Paper  I,  using  the  "FIVE” 
procedure.  The  proximity  on  the  sky  of  Saturn  and  the  reference  source, 
Jupiter,  permitted  one  improvement  over  the  previous  observing  procedures: 
measurements  of  Saturn  and  Jupiter  were  Interspersed  every  hour  or  so,  instead 
of  on  a  time  scale  of  many  hours  or  days.  This  fact  made  the  system  calibra¬ 
tion  and  atmospheric  attenuation  factors  for  the  measurements  of  the  two 
planets  even  more  uniform  than  previously,  thus  assuring  the  accuracy  of  the 
resulting  Saturn/ Jupiter  ratios.  Calibration  and  data  reduction  procedures 
were  also  as  before.  The  new  ratios  are  listed  and  displayed  in  Table  I  and 
Fig.  1,  respectively,  along  with  the  ratios  obtained  previously. 


» 


9 


RING  INCLINATION  ANGLE,  B 


Comparison  of  Che  observed  Saturn/Jupiter  ratios  with  the  simple  90- 
GHz  ring  model  of  Paper  I.  The  ratios  here  are  slightly  lower  than 
those  in  Table  I;  the  correction  for  the  size  of  Saturn's  disk  has 
been  removed  from  the  ratios  because  the  model  convolves  the  image  of 
Saturn  and  its  rings  with  the  Aerospace  antenna  pattern.  The  heavy 
curve  corresponds  to  the  simultaneous  least-squares  best-fit  to  the 
data  of  the  ring  brightness  temperature  T^  and  the  zero-inclination- 
angle  Intercept  of  the  Saturn/Jupiter  ratio.  The  thin  curves  indi¬ 
cate  the  1-sigma  limits  to  TR.  Utilizing  the  absolutely  calibrated 
Jupiter  brightness  temperature  Indicated  (Ullch  et  al.,  1980)  leads 
to  the  Saturn  disk  brightness  temperature  Tg  indicated. 


III.  ANALYSIS 


i 

As  in  Paper  I,  the  ring  brightness  temperature  has  been  derived  by  com¬ 
paring  the  data  with  the  variation  with  inclination  of  the  total  flux  from  the 
I  planet  and  rings  predicted  by  a  simple  model  with  uniformly  bright  A  and  B 

rings  (see  Figs.  3  and  4  of  Paper  I).  This  model  assumes  that  the  A  and  B 
ring  brightnesses  (TR)  are  equal,  uniform,  and  constant  with  inclination  angle 
B,  and  further,  that  the  C  ring  is  invisible.  The  disk  of  Saturn  is  assumed 
to  be  uniformly  bright  (Tg)  except  where  the  rings  cross  in  front  of  the  disk, 
in  which  case  the  brightness  is  the  sum  of  the  ring  brightness  plus  that  of 
the  disk  obscured  by  the  rings  (TR  +  Tg  exp[-Tg/sin(B)] ) ,  where  the  normal-to- 
the-ring  optical  depth  is  also  taken  to  be  uniform  across  the  rings  and 
equal  to  0.7.  The  free  parameters  of  this  model  are  the  brightness  tempera¬ 
tures  of  the  rings  (TR)  and  the  disk  (Tg),  and  are  determined  in  a  simul¬ 
taneous  least-squares  fit  of  the  data  to  the  model  variation  of  total  flux 
with  ring  inclination.  Taking  the  absolutely-calibrated  86-GHz  Jovian  disk 
temperature  and  its  1-sigma  uncertainty  Tj  »  177  ±  5  K  (Ulich  et  al.,  1980;  we 
multiplied  their  value  by  0.988  to  convert  it  to  the  Jovian  dimensions  used  in 
Paper  I  and  we  assumed  that  there  is  no  change  in  temperature  between  86  and 
90  GHz),  we  obtain 


The  errors  are  formal 
procedure,  and  include  the 
values  scale  directly  if  a 


Tr  -  17  ±  2.5  K 
Ts  -  156  ±  6  K. 

1-slgma  errors  from  the 
above  calibration  error 
different  flux  scale  is 


least  squares  fitting 
for  the  flux  scale.  Both 
chosen.  Figure  1  shows  a 
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comparison  of  the  data  with  the  model  result  for  the  best-fit  disk  and  ring 
temperatures  (heavy  curve).  The  light  curves  show  the  model  results  obtained 
using  the  1-sigma  limits  on  TR  and  the  best  fitting  Tg  for  each  case.  As 
predicted  in  Paper  1,  the  addition  of  the  low-inclination  points  greatly 
reduces  the  uncertainties  in  the  derived  values  of  TR  and  Tg;  e.g.,  the  uncer¬ 
tainty  in  TR  has  been  reduced  from  8  to  2.5  K,  while  Tg  is  essentially 
directly  measured  by  these  points. 

[Muhleman  and  Berge  (1984)  derived  TR  ■  13  ±  6  K  from  2.7-mm  inter¬ 
ferometry  in  October  1981,  which  agrees  with  our  result,  and  TR  *  32  ±  4  K 
from  interferometry  in  June  1983;  instrumental  problems  affect  both  results 
and  are  the  cause  of  the  discrepancy  (Muhleman,  1983).] 

While  this  model  is  certainly  oversimplified,  it  was  shown  in  Paper  I 
that  it  is  capable  of  considerable  generalization.  First,  the  model  is  insen¬ 
sitive  to  the  actual  distribution  of  brightness  within  the  rings.  By  solving 
for  TR  under  the  respective  assumptions  that  only  the  A  or  B  rings  contribute 
to  the  emission,  it  was  shown  that  the  product  of  TR  times  the  area  of  the 
contributing  ring  was  within  10  percent  of  that  for  the  case  of  uniformly 
bright  A  and  B  rings.  Hence  TR  can  be  considered  a  good  approximation  to  a 
mean  (area  averaged)  ring  brightness  temperature.  Secondly,  varying  the 
normal  optical  depth  from  0.4  to  1.0  resulted  in  a  total  variation  of  only  1.5 
K  in  TR,  so  that  the  result  is  insensitive  to  the  ring  optical  depth. 

The  effect  of  the  C  ring  was  neglected  in  Paper  I  because  of  its  small 
area  and  low  optical  depth.  We  reexamined  this  assumption  because  of  the  much 
smaller  error  on  our  present  result  for  the  ring  brightness  temperature,  and 
in  light  of  a  better  understanding  of  the  properties  of  the  C  ring  from 
Voyager  results.  Tyler  et  al.  (1983)  find  that  the  extinction  of  the  C  ring 
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is  about  0.2  at  3.6-cm  wavelength  from  1.22  Saturn  radii  to  the  inner  edge  of 
the  B  ring.  This  extinction  value  must  be  divided  by  2  to  obtain  the  optical 
depth  appropriate  for  the  present  radiative  transfer  problem,  because  strongly 
forward-scattered  radiation  still  contributes  to  the  brightness  of  the  disk 

seen  through  the  rings  (Cuzzi  and  Pollack,  1978,  and  Cuzzi  et  al.,  1980). 

Comparison  of  the  3.6-cm  optical  depth  with  that  at  visual  wavelengths  shows 
little  variation  beyond  this  effect  (Cuzzi  et  al.,  198A).  Assuming  that  the 
3-mm  optical  depth  is  the  same  as  that  at  3.6  cm,  we  included  the  C  ring  in 
our  model  with  a  normal  optical  depth  of  0.1  and  a  brightness  temperature 
modelled  as 

TRC>  "  eT0  (1  "  exPt“°* l/sintB| ] ) 

(C) 

where  TQ  is  the  physical  temperature  of  the  C  ring  particles  and  e  their 

average  emissivity.  Emission  from  the  region  of  the  disk  obscured  by  the  C 

ring  is  of  course  reduced  by  the  factor  expl-0. l/sin|B| ) . 

(C) 

The  C  ring  brightness  temperature  is  fully  determined  if  T^  and  e  are 
(C) 

chosen.  Taking  T^  '  ■  100  K  and  letting  e  take  the  extreme  values  of  0  and  1, 
we  again  solved  for  a  uniform  TR  for  the  combined  A  and  B  rings,  obtaining  TR 
-  19.2  and  13.9  K  respectively.  The  goodness  of  fit  by  a  chi-square  measure 
does  not  vary  significantly,  and  in  all  cases  is  fully  consistent  with  the 
measurement  errors.  If  we  consider  that  all  values  for  e  are  equally  likely, 
then  che  result  for  TR  originally  obtained  is  essentially  unchanged,  while  the 
uncertainty  Introduced  by  our  lack  of  knowledge  about  the  emissivity  of  the  C 
ring  is  comparable  to  that  of  our  original  formal  error  for  TR.  Conversely, 
the  present  results  tell  us  nothing  new  about  the  C  ring. 


We  conclude  Chat  Che  derived  ring  brightness  temperature  represents  a 
mean  brightness  temperature  for  the  A  and  B  rings  which  Is  valid  for  a  wide 
range  of  possible  models  for  the  entire  ring  system.  Taking  the  above  sources 
of  uncertainty  Into  account,  we  estimate  that  the  mean  ring  brightness  tem¬ 
perature  Is 

Tr  -  17  ±  4  K. 

This  result  presumes  that  the  ring  brightness  temperature  remains 
constant  as  the  ring  Inclination  varies.  However,  the  ring  brightness  will 
change  inasmuch  as  both  the  physical  temperature  of  the  ring  particles  and 
their  scattering  properties  will  depend  upon  inclination.  It  is  apparent  from 
Fig.  4  of  Paper  I  chat  the  model  is  most  sensitive  to  TR  for  B  £  10  deg  and  is 
not  sensitive  to  TR  at  low  Inclinations  (where  the  largest  brightness  varia¬ 
tions  would  be  expected  to  occur)  so  that  the  present  result  is  weighted 
toward  the  larger  inclinations.  A  more  detailed  model  is  needed  to  make  this 
consideration  more  quantitative  -  such  a  model  is  described  in  Sec.  IV.  B 
below  and  leads  to  a  ring  brightness  temperature  averaged  over  radius  and 
inclination  which  is  consistent  with  the  above  result. 

iv.  utekpretation 


We  use  the  same  two  methods  outlined  in  Paper  I  to  obtain  improved  limits 
on  particle  and  ring  properties  on  the  assumption  that  the  particles  are  pure 
water  ice. 
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L.  "Thin  Slice  of  9irtti  Cheese*  Model 

Diffuse  scattering  of  the  planet's  thermal  emission  produces  a  radially 
dependent  component  of  the  total  ring  brightness  of  average  magnitude  6  ±  3  K 
(Paper  I).  Subtracting  out  this  scattered  emission  leaves 


Tt  -  11  ±  5  K. 


We  regard  this  latter  component  as  representing  the  intrinsic  thermal  emission 
from  the  rings;  it  is  now  a  definite  result  rather  than  the  upper  limit 
obtained  in  Paper  I.  This  result  allows  us  to  estimate  the  mean  rather  than 
the  limiting  properties  of  the  rings •  For  example,  consider  the  simple  slab 
(or  “thin  slice  of  Swiss  cheese")  model  (Paper  I).  For  this  model  the  rings 
are  regarded  simply  as  a  uniform  rough  slab  (but  with  holes  in  it  to  represent 
the  spaces  between  particles)  of  thickness  r,  normal-to-the-ring,  optical 
depth  Tq,  and  physical  temperature  Tq;  the  thermal  component  Tt  of  the  bright¬ 
ness  temperature  is  given  approximately  by 


T 


t 


oar  (l 


(1) 


where  a  is  the  wavelength-dependent  absorption  coefficient  of  the  ring 
material.  The  exponential  term  accounts  for  the  open  area  between  particles 
and  is  a  good  approximation  for  a  monolayer  if  B  is  not  too  small  (Cuzzie  and 
Pollack,  1978).  We  take,  as  in  Paper  I,  Tq  csc|b|  *0.7  esc  20°  *  2  and 
Tq  ■  85  K  (derived  from  22. 7- ps  ground-based  observations  at  B*  *  16  and  26* 
(see  No It  et  al. .  1978)].  From  extensive  new  laboratory  measurements  of  pure 
water  ice,  Mishima,  Klug,  and  Whalley  (1983)  have  derived  at  3.3  mm  and  85  K, 
an  absorption  coefficient  a  -  (4.8  t  1.0)  x  10  Jcm  .  (This  value  is  1.5 
times  larger  than  the  value  adopted  In  Paper  I.)  We  then  obtain 
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r  *  30  ±  13  cm 


for  Che  Chickness  of  a  simple  ice-slab  model  ring.  Because  Tq  decreases 
somevhaC  as  |B|  decreases,  this  value  of  r  is  possibly  somewhat  low.  A  more 
detailed  model  is  discussed  next. 

B.  Multiple-Scattering  Model 

This  method  employs  a  numerical  radiative  transfer  model  in  which  the  A, 
B,  and  C  rings  are  represented  by  scattering/emitting  layers  of  finite  optical 
depth.  The  exact  particle  single-scattering  albedo  (u^>  and  phase  function 
are  determined  by  Mie  scattering  calculations,  assuming  only  the  particle  size 
distribution  and  refractive  indices  nr  +  in^  -  »r  Mfe  a]l.  Because  of  the 
ensuing  high  particle  albedos,  the  net  brightness  of  the  rings  results  from  a 
combination  of  ring  thermal  emission  and  diffuse  reflection  of  planetary 
thermal  emission.  This  model  is  described  in  Paper  I  and  in  more  detail  by 
Cuzzl  et  al.  (1980);  it  is  compared  with  the  observations  in  Fig.  2.  The 
heavy  curve  corresponds  to  the  simultaneous  least-squares  best  fit  to  the  data 
of  the  zero-inclination  Intercept  of  the  Saturn  system/Jupiter  ratio 
and  lsot,  which  is  a  scaled  effective  albedo  for  isotropic  scattering 
obtained  [by  similarity  relations  (Paper  I  and  Irvine,  1975)]  from  the  real 
albedo  and  a  (highly  forward-scattering)  phase  function.  The  thin  curves 
correspond  to  the  1 -sigma  limits  on  iBOt*  best-fit  values 

°e  *b,i.ot  *"d  Ts  *re 


ViMt  ■  °-85  *  °-°5 


Tg  -  156  *  6  K, 
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RING  INCLINATION  ANGLE,  B 


Fig.  2  Comparison  of  the  observed  90-GHz  Saturn/Jupiter  ratios  with  the 

multiple-scattering  model.  The  heavy  curve  corresponds  to  the  simul¬ 
taneous  least-squares  best-fit  to  the  data  of  the  equivalent  iso¬ 
tropic  single-scattering  albedo  and  the  aero-inclination- 

angle  intercept  of  the  Saturn/Jupiter  ratio.  The  thin  curves  indi¬ 
cate  the  1-  sigma  limits  to  u>q  igot*  Utilizing  the  absolutely  cali¬ 
brated  Jupiter  disk  brightness  shown  (Ulich  et  al.,  1980)  leads  to 
the  Saturn  disk  brightness  temperature  Tg  indicated. 


where  the  error  for  Tg  accounts  for  uncertainties  in  the  ring  optical  depth 
(which  are  negligible — see  Paper  I  and  Sec.  Ill),  in  Jupiter’s  disk  tempera¬ 
ture,  and  in  the  fitting.  In  these  calculations  the  ring  physical  temperature 
is  assumed  to  vary  with  solar  elevation  angle  B',  as  is  observed  [see  the 
discussion  in  Paper  I  and,  e.g.,  Nolt  et  al.  (1978)].  The  inclination- 
averaged  ring  brightness  temperature  computed  using  this  multiple-scattering 
model  and  igot  ■  0.85  is  *  15  K  —  consistent  with  the  value  of  17  ±  3  K 
derived  from  the  simpler  model  of  Sec.  Ill  which  assumes  the  ring  brightness 
to  be  Independent  of  Inclination  angle. 

We  also  investigated  the  possible  effect  on  the  results  of  the  presence 
of  "holes"  in  the  rings.  If  a  few  percent  of  the  ring  area  has  an  optical 
depth  of  only  a  few  percent  [as  implied  by  analyses  of  Pioneer  Saturn  data 
(Esposito  et  al.,  1980)],  then  blockage  of  the  planetary  emission  will  be  less 
at  low  inclination  angles  than  if  a  homogeneous  ring  were  assumed.  Such 
"holes"  will  cause  our  multiple  scattering  model  to  derive  a  spuriously  low 
value  of  lBOt«  but  by  an  amount  less  than  the  above  formal  error.  Thus  if 
"holes"  exist,  the  true  value  of  may  be  slightly  higher  than  quoted 

above.  Voyager  observations,  however,  do  not  support  such  a  large  fraction  of 
empty  space  (Lane  et  al.,  1982;  Esposito  est_al_. ,  1983).  Consequently,  in  this 
regard,  our  value  of  igQt  is  acceptable. 

The  variations  of  the  particle  single-scattering  albedo  u>Q  and  the  iso¬ 
tropic  single-scattering  albedo  ^got  as  a  function  of  absorption  coeffi¬ 
cient  o,  assumed  size  distributions  for  the  particle  volume  density  n(r)  ■ 
nQr"8,  and  maximum  particle  size  rMX  (assuming  raln  •  1.0  cm  for  all  cases) 
are  shown  in  Table  II  and  Fig.  3.  The  tabulated  values  of  and  the 
asymmetry  parameter  g  are  obtained  by  averaging  over  the  Indicated  particle 
distribution;  lgQt  is  obtained  from  and  g  (see  Table  II,  Footnote  c). 
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TABLE  II 


VARIATION  OF  THE  PARTICLE  SINGLE-SCATTERING  ALBEDO  ZQ  AND  THE  ISOTROPIC 

SINGLE-SCATTERING  ALBEDO  .  OF  PURE  WATER-ICE  Ih  PARTICLES  WITH  SIZE 

u,xsot 

DISTRIBUTION  n(r)  »  rigr"8  FROM  r  -  1  cm  TO  r  -  r mMX  AS  A  FUNCTION  OF  s, 

AND  ABSORPTION  COEFFICIENT  a/a  8 

non 


• 

r«M 

(a) 

a/ a  * 

no* 

; b 
“o 

gb.C 

-v  b,C 

Li 

0»l«ot 

3.0 

1.0 

1.0 

0.954 

0.749 

0.838 

3.0 

5.0 

1.0 

0.869 

0.775 

0.599 

3.0 

1.0 

0.3 

0.963 

0.738 

0.940 

3.0 

3.0 

0.3 

0.953 

0.749 

0.840 

3.0 

5.0 

0.3 

0.926 

0.758 

0.750 

3.0 

1.0 

0.1 

0.994 

0.734 

0.980 

3.0 

3.0 

0.1 

0.981 

0.739 

0.930 

3.0 

5.0 

0.1 

0.967 

0.744 

0.880 

3.3 

1.0 

1.0 

0.968 

0.737 

0.888 

3.3 

3.0 

1.0 

0.931 

0.7*7 

0.780 

3.3 

5.0 

1.0 

0.906 

0.755 

0.700 

3.3 

1.0 

0.3 

0.988 

0.729 

0.960 

3.3 

3.0 

0.3 

0.968 

0.736 

0.890 

3.3 

5.0 

0.3 

0.949 

0.7*3 

0.830 

3.3 

1.0 

0.1 

0.996 

0.726 

0.990 

3.3 

3.0 

0.1 

0.987 

0.730 

0.950 

3.3 

5.0 

0.1 

0.978 

0.733 

0.920 

3.6 

1.0 

1.0 

0.977 

0.727 

0.921 

3.6 

3.0 

1.0 

0.954 

0.733 

0.850 

3.6 

5.0 

1.0 

0.937 

0.738 

0.796 

“noe 

-  6.6  «  IQ"3 

c*~*  at  3.3  m  (Mlahlaa  at  al. 

.  1983). 

b  Theta  valuta  art 

dltCribatlon. 

obtained  by  averaging 

over  the 

Indicated  particle  elte 

c  Th« 

laotroplc  alngle-ecatterlng  albedo 

*t,laot 

l0u  -  g)/U 

-  “o8>- 

where  the  aaywwetry  ptriMttt  g  •  2  «  cm 9  p(8)  ein#  d6  (aee  Paper  I). 
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MAXIMUM  PARTICLE  SIZE,  rmax  Iml 

Fig.  3  Relationships  between  the  isotropic  single-scattering  albedo  igQt  at 
3.3  no  and  the  value  of  rMX  for  an  ensemble  of  pure  water-ice  par¬ 
ticles  with  size  distribution  n(r)  -  i:)r“8 ,  (s  ■  3.0,  3.3,  and  3.6) 
from  r  ■  1  cm  to  r  *  tMX  and  with  a  3.3-mm  absorption  coefficient 

equal  to  that  of  pure  water  ice  Ih  at  85  K  (a/a  -  1.0).  Also  shown 

non 

are  the  relationships  for  particles  with  absorption  coefficients  0.1 
and  0.3  that  of  pure  water  ice  ( a/ a  ■  0.1  and  0.3).  The  best-fit 

nOm 

value  of  “  0.85  ±  0.05,  as  derived  from  comparison  of  the  3.3- 

mm  observations  with  the  multiple-scattering  model,  is  Indicated  by  the 
cross-hatching. 
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The  observed  absorption  coefficient  a  is  normalized  by  the  nominal  value  at 

3.3  mm  for  pure  water  ice  Ih  at  85  K  of  a  ■  4.8  x  10”^  cm"*  (Mishima  et 

nom  — 

al. ,  1983).  Values  used  for  the  slope  s  in  the  power  law  are  3.0,  3.3,  and 

3.6;  this  range  is  consistent  with  ground-based  radar  observations  (Cuzzi  and 

Pollack,  1978).  Further,  Voyager  results  indicate  that  the  slope  s  is  between 

3.2  and  3.4  (Marouf  et  al..  1982).  Marouf  et  al.  also  deduced  that 

rmax  ~  5  m  in  most  of  the  ring  system;  only  in  the  C  ring,  which  contributes 

little  to  the  3-mm  emission,  is  rmax  possibly  as  small  as  1  m.  Inspection  of 

Fig.  3,  where  our  result  for  uL  .  _  of  0.85  ±  0.05  is  indicated  by  the  cross- 

U | X80t 

hatching,  shows  that  for  our  result  to  be  consistent  with  an  rmax  as  large  as 
5  m,  either  the  slope  s  must  be  ~3.6,  or  the  particle  absorption  coeffi¬ 
cient  a  at  3  mm  must  be  0. 1  to  0.3  of  a  ,  or  some  combination  of  the  two.  A 

nom 

slope  of  3.6  is  somewhat  larger  than  the  best-fit  Voyager  result,  but  not 

Impossibly  so.  Since  the  ring  particles  are  unlikely  to  be  pure  water  ice 

(Clark,  1980),  or  since  they  could  contain  ice  in  its  amorphous  phase,  it  is 

possible  that  the  particles'  absorption  coefficient  could  be  lower  than  that 

of  pure  ice  Ih.  The  calculations  displayed  in  Fig.  3  indicate  that  self- 

consistency  between  all  the  constraints  is  achieved  if  rmax  ~  5  m,  s  ~  3.3, 

and  a  ~  0.5  a 

nom 

C.  Caveat 

If  lnterpartlcle  shadowing  causes  particles  deep  within  the  ring  layer  to 
be  physically  cooler  than  those  which  are  unshadowed,  then  some  of  this  dis¬ 
crepancy  In  a  may  be  removed.  Furthermore,  we  note  that  the  Voyager  size 
distribution  refers  to  regions  of  low  to  moderate  optical  depth,  whereas  our 
measurements  refer  to  regions  of  moderate  to  high  optical  depth  where  the 
exact  size  distribution  is  not  as  well  known  and  may  be  different. 


22 


We  also  Investigated  the  effect  of  merely  decreasing  the  density  of  the 


ring  particles.  The  refractive  index  was  assumed  to  vary  with  density  as 
predicted  by  the  Rayleigh  mixing  formula,  as  suggested  by  laboratory  measure¬ 
ments  (Campbell  and  Ulrichs,  1969).  The  effect  of  decreasing  density  is  thus 
to  decrease  both  the  real  and  imaginary  indices  of  refraction.  The  main 
effect  on  particle  scattering  properties  was  to  increase  the  degree  of  forward 
scattering,  but  not  to  change  the  albedo.  Model  calculations  of  ring  tempera¬ 
tures  showed  that  the  fit  to  the  brightness  temperature  data  was  not  Improved 
by  such  a  decrease  in  particle  density,  and  that  implications  for  rmax  and  s 
were  not  affected.  Therefore,  we  feel  that  there  is  no  evidence  that  the 
particle  density  is  much  lower  than  that  of  solid  ice.  However,  these  data 
cannot  rule  out  the  presence  of  a  particle  density  somewhat  less  than  unity. 

V.  DISCUSSION 


The  low  brightness  temperature  strongly  implies  the  presence  of  a  low- 
loss  dielectric  and  places  a  limit  on  the  amount  of  high-loss  material  which 
can  be  present.  In  the  present  case  we  may  estimate  an  upper  limit  to  the 
rock/lce  ratio  within  the  rings.  The  absorption  coefficient  of  chondritic 
type  material  at  3.3  mm  ranges  upward  from  approximately  0.8  cm-1  (Campbell 
and  Ulrichs,  1969).  At  lower  densities  (£»£.,  loose  powders  or  dispersed 
particles),  the  absorption  coefficient  of  a  mixture  is  approximately  propor¬ 
tional  to  the  fraction  of  the  absorber.  Hence,  from  Eq.  (1),  distributed 
material  of  total  quantity  equivalent  to  a  slab  of  2  mm  thickness  would 
produce  a  3.3-mm  brightness  tempersture  of  15  K  for  Tq  -  85  K.  This  Is  the 
largest  quantity  of  such  rocky  material  which  can  be  present,  if  uniformly  and 
finely  divided.  On  the  other  hand,  the  minimum  quantity  of  Ice  required  in 
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order  to  scatter  radar  at  X-  and  S-band  is  a  monolayer  of  particles  of 
radius  >  6  cm  (Cuzzi  and  Pollack,  1978).  Therefore,  allowing  for  a  rock 

density  of  2.4  g/cm^,  the  rock/ice  ratio  must  be  less  than  about  10%  if  the  ( 

materials  are  uniformly  mixed,  consistent  with  radar  observations  (Cuzzi  and 

Pollack,  1978).  A  more  stringent  limit  may  be  placed  using  our  derived  values 

of  ~  5  m,  s  ~  3.3  and  a  ~  0.5a  (i.e.,  the  observed  absorption  coeffi-  . 

max  nom  1 

cient  a  is  about  half  the  nominal  absorption  coefficient  of  ice).  If  the 
rocky  material  were  a  minor,  uniformly  distributed  fraction  and  the  dominant 

icy  material  had  an  absorption  coefficient  of  zero,  there  could  be  no  more  ( 

than  (observed  a)/(rock  a)  ~  2.4  x  10“^/0.8  ~  3  x  10“^  by  volume  of  any  rocky 
or  carbonaceous  absorber. 

This  result  argues  strongly  for  some  sort  of  chemical  differentiation  ( 

having  taken  place  within  the  rings.  Several  current  theories  for  the  origin 
of  the  rings  involve  condensation  more  or  less  in  situ  from  a  cloud  of  solar 
composition  around  Saturn  (e.£. ,  Pollack,  1975;  Cameron,  1977).  Without  dif¬ 
ferentiation,  and  assuming  that  all  the  H2O  (but  not  necessarily  all  the  NH3 
or  CH^)  is  retained,  then  the  ratio  of  rock/ ice  would  be  about  50% 

(Consolmagno  and  Lewis,  1977).  The  lowest-density  Saturnian  satellites  seem 
to  have  about  a  20%  rock/lce  ratio  (Smith  et  al«,  1982).  A  ratio  of  10%  or 
less  would  seem  to  imply  that  an  even  more  significant  loss  of  rock  has  taken 
place  In  the  rings  than  in  the  satellites  relative  to  the  original  solar 
nebula. 

However,  it  is  not  necessary  to  physically  remove  the  rock  completely 
from  the  rings  to  achieve  the  present  result.  Rather,  the  rock  need  only  be 
hidden  from  view  at  microwave  frequencies.  For  example,  we  may  consider  a 
blmodal  distribution  of  material  within  the  rings  such  that  the  rock  is  con¬ 
tained  in  a  small  number  of  relatively  large  particles,  while  the  ice  is  more 
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broadly  distributed,  possibly  in  different  particles.  If  the  rock  particles 
are  a  few  meters  or  larger  in  diameter,  then  they  may  contain  an  effectively 
unlimited  mass  relative  to  the  ice  without  contributing  significantly  to  the 
thermal  emission  because  their  net  cross-sectional  area  remains  small. 

Such  a  situation  could  arise  as  follows:  Differential  condensation  of  a 
I  solar  composition  cloud  around  Saturn  could  have  occurred  according  to  the 

model  of  Consolmagno  and  Lewis  (1977),  leaving  initially  a  swarm  of  small 
rocky  bodies  of  meter-size  and  larger.  Subsequent  cooling  of  the  nebula  would 
|  result  in  further  condensation  of  volatiles,  producing  an  ice  coating  on  these 

bodies.  As  they  are  within  Roche's  limit  (i.e.,  in  the  vicinity  of  the 
rings),  these  "moonlets"  may  not  have  accreted  further  into  larger  bodies,  but 

► 

may  remain  more  or  less  in  place  restricted  in  size  («  1  km)  at  the  pre¬ 
sent.  Meanwhile,  repeated  weak  collisions  among  these  bodies  have  chipped  the 
ice  from  their  surfaces.  These  ice  particles,  with  sizes  of  centimeters  to 
meters,  produce  most  of  the  visible  surface  of  the  rings.  The  surface  area  of 
the  parent  objects  is  much  less  than  the  surface  area  of  the  ’’visible"  ring 
particles.  Furthermore,  their  total  surface  area  is  far  less  than  if 
uniformly  dispersed  Into  small  particles.  This  model  requires  only  that 
F  "local"  differentiation  has  taken  place.  Alternatively,  of  course,  dlfferen- 

(  tlatlon  processes  on  larger  scales  may  have  been  active  to  segregate  the  rocky 

^  material  or  even  to  remove  this  material  entirely  from  the  vicinity  of  the 

I  rings,  e.g.,  through  viscous  drag  on  the  initial  condensation  of  rocky 

| 

moonlets  by  the  as-yet-uncondensed  volatile  component  (Pollack,  1975).  The 
^  anomalously  low  densities  of  Saturn's  satellites  measured  by  Voyager  could 

I  also  be  Indicative  of  such  a  process.  Preliminary  photometry  of  Voyager 

imaging  results  points  to  variations  of  particle  albedo  in  the  rings 

k 

I  from  ~  0.3  in  the  C  ring  to  ~  0.6  In  the  A  and  B  rings;  also,  closer 
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inspection  shows  comparable  variations  in  albedo  even  on  scales  of  50  -  100  km 
in  the  B  ring  (Cuzzi  et  al.,  1984).  This  inhomogeneity  is  supportive  of  our 
suggestions  above. 
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LABORATORY  OPE RAT I ONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting  exper¬ 
imental  and  theoretical  Investigations  necessary  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  military  space  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  in  the  nation's  rapidly  developing 
space  systems.  Expertise  In  the  latest  scientific  developments  Is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are; 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high-temperat ure  thermoraechanlcs ,  gas  kinetics  and  radiation;  research 
in  environmental  chemistry  and  contamination;  cw  and  pulsed  chemical  laser 
development  including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bioenvironment al  research  and 
monitoring. 

Electronics  Research  Laboratory:  Hie  coelectronics,  GsAs  lov-nolae  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  lidar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  mi lliaeter-wave  and  microwave  technology. 

Information  Sciences  Research  Office:  program  verification,  program  trans¬ 
lation,  performance-sensitive  system  design,  distributed  architectures  for 
apaceborne  computers,  fault-tolerant  computer  systems,  artificial  intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analvsls  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
space  environment;  materials  performance  In  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  In  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  phvsics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  on  apace  systems. 


